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Abstract 

It is shown that an event sample from the Monte Carlo simula- 
tion of a random cascading a model with fixed dynamical fluctuation 
strength is intermittent but not chaotic, while the variance of dy- 
namical fluctuation strength in different events will result in both the 
intermittency and the chaoticity behavior. This shows that fractality 
and chaoticity are two connected but different features of non-linear 
dynamics in high energy collisions. 



^This work is supported in part by the National Natural Science Foundation of China. (NSFC) 
under Grant No.l9575021. 



For classical system, the description of non-linear behavior is well established. It has 
been known by lattice calculation that the classical non-Abelian gauge theory gener- 
ally exhibits deterministic chaos and that the Lyapunov exponent can be numerically 
determined But for quantum system, because of the ambiguity associated with 
quantum chaos in the realm of quantization, nonconservation of the number of degrees 
of freedom and lack of a meaningful definition of a trajectory, there is no correspond- 
ing description existed. The study of non-linear behavior in high energy physics has, 
therefore, to be started phenomenologically. 

The first signal of such a behaviour came from the unexpectedly large local fluctua- 
tions in a single event of very high multiplicity recorded by the JACEE collaboration 0] . 
Such large fluctuations may not be simply due to statistical reason and was taken as a 
signal of the existence of non-linear dynamical fluctuations. It was soon realized that 
the idea can be applied to events of any multiplicity provided that a proper averag- 
ing of factorial moments is performed, as done in the pioneer work of Bialas and 
Peschanski. These authors have been able to show that, if the statistical fluctuations 
are of Bernouli (fixed multiplicity case) or Poisson (variable multiplicity case) type, 
the averaged factorial moments Fg is equal to the averaged dynamical probability mo- 
ments Cq. The anomalous scaling of the latter has taken the name of intermittency (or 
fractal). This led to extensive experimental studies and the expected anomalous 
scaling has been observed successfully in the experiments 0. 

It should be realized, however, that the averaging procedure, apart from its clear 
advantages, brings also a danger of losing some important information on the spatial 
patterns from event to event. In particular, some interesting effects, if present only in 
a part of events produced in high-energy collisions, may be missed. A possible example 
of this kind is the quark-gluon plasma which is expected to be characterized by specific 
intermittency exponents 0. It seems therefore important and urgent to study the 
fluctuation of single-event moments C'^'^ inside an event sample^. This fluctuation is 
related to the chaotic behavior of the system |]^. A quantity iiq called entropy index 
can be introduced p as an adequate parameter in measuring the chaotic behavior. 
The positivity of entropy index /i^ > is proved to be a criterion for chaos @. 

Thus, two kind of non-linear phenomena — fractality (intermittency) and chaoticity 
have been proposed in high energy collisions. In this short note we will study the 
relation beteen them using Monte Carlo simulation of a random cascading a model. We 
will show that the anomalous scaling of the averaged probability moments (fractality or 
intermittency) and that of the event-space moments of single-event ones (chaoticity) are 
two connected but different features of non-linear dynamics. The system will exhibite 
both the fractal and the chaos behaviour only when the dynamical fluctuation strength 
is not fixed but is distributed over a certain range. 

Let us first recall briefly the study of fractality (intermittency) in high energy 
collisions. This study is performed through the observation of anomalous scaling of 
averaged factorial moments F^, which is equal to the averaged probabilty moments Cq 

F,(M) = C,(M)^-E^«M^^ (1) 

^It has been shown [|| that the statistical fluctuations can be ehminated by using the factorial 
moments averaged over event sample. However, the extension of this method to single-event moments 
is highly non-trivial. It is easy to show that the elimination of statistical fluctuations in single- 
event factorial moments Fq is incomplete. In order to avoid the complication caused by statistical 
fluctuations we will in this paper restrict ourself to the study of probability moments C^"^ directly. 



where a phase space region A is divided into M sub-cells, pi is the probability for a 
particle to fall in the ith sub-cell. 

For a flat inclusive distribution the moment C'f'^ for each event is defined as 
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We can now consider C^^-* not only through its average — intended to get a better 
estimate of the hypothetical anomalous scaling of single-bin moments, cf. eqn.(l) — 
but also as a pattern-descriptor for particle fluctuations inside bins (just one among 
the many that could be devised). 

C'^'^ may fluctuate greatly from event to event. In a sample consisting of a large 

number of events, we get a distribution of Cf>, denoted by P {C^f^^, which is 
normalized to unity. The conventionally defined factorial moments, cf eqn.(l), give 
only an estimate of the mean of P {C^^^. By taking the normalized moments of 

P (cj^^^^ in event-space defined as 

^p,. = (c^^/(c^^^)^ (3) 

we have a quantification of the fluctuation of the spatial patterns, i.e. we can investigate 
the full shape of the distribution and, especially, the way it changes with the resolution 
6 = A/M. The value of p can be any positive real number. If Cp^q{M) has a power 
law behaviour in M, i.e. 

Cp,g(M) oc M^^^P\ (4) 
then a new entropy index can be defined as. 



(5) 

p=i 



It is easy to see that finite, nonvanishing positive values of fig corresponds to wide 

which in turn means unpredictable spatial pattern from event to event. By 

applying the measure to known classical chaotic system, it has been shown that fig 
can be used as a measure of chaos in problems where only the spatial patterns can be 
observed and the positivity of fig is a criterion for chaos. 

An alternative way of calculating fig is to express Cp^g as 

Cr,,= {^f^'), (6) 

in which, 

$(^) = Cf /(Cf). (7) 

With the definition 

S, = (<|.f)ln<|.(^)), (8) 



we can obtain 



in the scaling region, i.e. where exhibits a hnear dependence on InM. We will use 
this formula in calculating the entropy indices Hg. 

Let us turn now to the consideration of the relation between the fractality and 
chaoticity in high energy collisions. Since the random cascading a-model 0|TD[ is 
often used to study the dynamical fluctuations in these collisions, we will use this 
simple model tool for our investigation. 

In the random cascading a-model, the M divisions of a phase space region A are 
made in steps. At the first step, it is divided into two equal parts; at the second step, 
each part in the first step is further divided into two equal parts, and so on. The steps 
are repeated until M = AY/6y = 2^ . How particles are distributed from step-to-step 
between the two parts of a given phase space cell is defined by independent random 
variable ujyj^^ where is the position of the window (1 < jy < 2'^) and u is the number 
of steps. It is given by UTOf : 



^u,2j-i = ^(1 + «r) ; ujy^2j = ^{l-ar), j = 1,...,2'' ^ 

where, r is a random number distributed uniformly in the interval [—1,1]. a is a 
positive number less than unity, which determines the region of the random variable 
uj and describes the strength of dynamical fluctuations in the model. After z/ steps, 
the probability in the mth window (m = 1, . . . , M) is Pm = '^iji^2j2 • • • ^vj^- Then 
according to eqn(l), probabilty moment C^*^-' in each event of different division steps 
are calculated, and the moment Cp^g and entropy index fig of the sample are obtained 
using eqn.(3) and eqn.(9). 

Our research is done in the following two steps. 

(A) Fix the model parameter to a definit value, say a = 0.34, the experimental 
results being around this value. The results of InCg, In Cp^g and Eg vs. InM from 6000 
MC simulation events are shown in Fig. 1(a), (6), (c) respectively. 

In Fig. 1(a) we see a straight line in bi- logarithm plot, which is an indication of 
intermittency (or fractality). However, the behavior of In Cp_g vs. InM in the model, 
cf. Fig. 1(6), is much different from the expected result for chaos [0]. It does not 
show any scaling behavior or upward bending when M goes larger as the chaotic 
behaviour requires [0 . The first going up of Cp^g is due to an intrinsic uncertainty of the 



intermittency parameters |]11|. The cascade responsible for intermittent behaviour has 
different realizations in different events, and the intermittency exponents determined 
from different realizations of the same random cascade are scattered around the average, 
i.e. the method has a finite resolution with respect to the parameters of the random 
cascade. The Cp^g saturates when M goes large means that there isn't any essential 
fluctuation of spatial pattern from event to event. Therefore, this kind of a model 
cannot reflect the feature of chaoticity. From the result showing in Fig. 1(c), using 
eqn.(9), we can get for this case fig ~ 0. If we donot consider the finite resolution of 
model, there isn't any chaotic behavior. 

In order to reproduce both intermittency (fractal) and chaos we take the second 
step. 

(B) Instead of giving a a fixed value we let it be a random variable having a 
Gaussian distribution. The mean value and variance of the Gaussian are both chosen 
as 0.22. Calculating from 6306 events, the result of InCg vs. InM are shown in Fig.2. 
It can be seen from the figure that there is very good power-law or scaling behavior, 
which means that though we have changed the method of setting model parameter, 
the anomalous scaling of the mean value of C^^-* (intermittency phenomenon) survives. 



With the method developed in Ref. [|T2| we can get the effective ffuctuation strength 
in this case as = 0.337. This value of a^s is within the limited range available in 
actual experiments. However, the behavior of In Cp^q vs. In M in the present case is 
much different from the case (A) and shows a typical behaviour of chaoticity, cf. Fig.3. 

From this result we see that a distribution of a will cause a distribution of single 
event probabilty moment Cj^'^\ i.e. for an event sample we have a wide P{C'^^). For 

increasing M, along with the expected increase of the average, P{Cf>) will show a 

rapid broadening, i.e. a more violent fluctuation of C'^^M) for different events, which 
will result in a even more unpredictable spatial pattern from event to event. 

Using eqn.(8) and eqn.(9), we can calculate entropy indices for this case. The result 
of Tjq vs. InM are shown in Fig.4(a). By performing a linear fit of vs. InM in 
the range M = 8 to M = 64 (i.e. omitting the first three points), /ig is obtained and 
plotted in Fig. 4(6). 

By this two steps of MC simulation we can see that the procedure of doing simula- 
tion with random cascading a model of fixed strength parameter a, as has been widely 
used before, captured only one aspect of the non-linear property (intermittency) but 
cannot reproduce the fluctuation of spatial patterns from event to event. It will cause 
the losing of information on the spatial patterns in different events and some intersting 
effects if they are present only in a part of events, may be lost. If we want to give 
a more complete description of the non-linear properties using a model, the model 
parameter cannot be fixed. 

In conclusion, the nonvanishing positive values of /Xg, which is an indication of 
chaos, correspond to wide P(C^'^^), which in turn means unpredictable spatial pattern 

from event to event. Such an unpredictability of wide P{C'^'>) is caused by different 
dynamical fluctuation strength in different events, i.e. by a distribution of dynamical 
fluctuation strength in an event sample. Events in one sample are all beginning with 
similar initial condition. During the collision process each event will evolve with a 
different strength of dynamical fluctuation and result in a fluctuation of spatial pattern 
in final state event space. 

Dynamical fluctuation strength is directly related to the dynamical mechanism in 
a particular collision. We take the distribution of dynamical fluctuation strength, i.e. 
the distribution of model parameter a, to be a Gaussian only because it is the most 
common distribution of random variables in nature. We have also tried a uniform 
distribution of a and non-zero positive entropy indices yUg can be obtained too. (The 
result is not shown here). Revealing the distribution of dynamical fluctuation strength 
in different collisions will be a very constructive work and it will certainly help us a lot 
in studying the mechanism of strong interactions. How the different distributions of 
dynamical fluctuation strength together with it's mean value and width will influence 
the entropy index of an event sample is also a problem worthwhile further investigation. 

As has been stressed in the introduction, this study is restricted to the probability 
moments and the problem of eliminating statistical fluctuations in experimental data 
analysis has been postponed. To develope an effective method for eliminating the 
statistical fluctuations for single-event moments is a challenge for future investigation. 
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Figure Captions 

Fig.l (a) Averaged C2, (b) Cp.2, (c) E2 for fixed a. Full lines are for guiding the eye. 

Fig.2 Averaged C2 for Gaussian-distributed a. 

Fig.3 The InCp^^ vs. InM for the random cascading model with 
Gaussian-distributed a. Full lines are for guiding the eye. 

Fig.4 Eg and /iq for Gaussian-distributed a. Full hues are for guiding the eye. 




Fig.l (a) Averaged C2, (6) Cp,2, (c) S2 for fixed a. 
Full lines are for guiding the eye. 
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Fig. 2 Averaged C2 for Gaussian-distributed 




Fig. 3 The InCp^g vs. InM for the random cascading model with 
Gaussian-distributed a. Full lines are for guiding the eye. 




Fig. 4 Tiq and iJ,q for Gaussian-distributed a. 
Full lines are for guiding the eye. 



